home *** CD-ROM | disk | FTP | other *** search
/ io Programmo 60 / IOPROG_60.ISO / soft / c++ / gsl-1.1.1-setup.exe / {app} / src / cblas / test_copy.c < prev    next >
Encoding:
C/C++ Source or Header  |  2002-04-18  |  4.7 KB  |  222 lines
  1. #include <gsl/gsl_test.h>
  2. #include <gsl/gsl_ieee_utils.h>
  3. #include <gsl/gsl_math.h>
  4. #include <gsl/gsl_cblas.h>
  5.  
  6. #include "tests.h"
  7.  
  8. void
  9. test_copy (void) {
  10. const double flteps = 1e-4, dbleps = 1e-6;
  11.   {
  12.    int N = 1;
  13.    float X[] = { 0.898f };
  14.    int incX = 1;
  15.    float Y[] = { 0.699f };
  16.    int incY = -1;
  17.    float expected[] = { 0.898f };
  18.    cblas_scopy(N, X, incX, Y, incY);
  19.    {
  20.      int i;
  21.      for (i = 0; i < 1; i++) {
  22.        gsl_test_rel(Y[i], expected[i], flteps, "scopy(case 76)");
  23.      }
  24.    };
  25.   };
  26.  
  27.  
  28.   {
  29.    int N = 1;
  30.    double X[] = { 0.002 };
  31.    int incX = 1;
  32.    double Y[] = { -0.921 };
  33.    int incY = -1;
  34.    double expected[] = { 0.002 };
  35.    cblas_dcopy(N, X, incX, Y, incY);
  36.    {
  37.      int i;
  38.      for (i = 0; i < 1; i++) {
  39.        gsl_test_rel(Y[i], expected[i], dbleps, "dcopy(case 77)");
  40.      }
  41.    };
  42.   };
  43.  
  44.  
  45.   {
  46.    int N = 1;
  47.    float X[] = { -0.166f, 0.639f };
  48.    int incX = 1;
  49.    float Y[] = { 0.863f, 0.613f };
  50.    int incY = -1;
  51.    float expected[] = { -0.166f, 0.639f };
  52.    cblas_ccopy(N, X, incX, Y, incY);
  53.    {
  54.      int i;
  55.      for (i = 0; i < 1; i++) {
  56.        gsl_test_rel(Y[2*i], expected[2*i], flteps, "ccopy(case 78) real");
  57.        gsl_test_rel(Y[2*i+1], expected[2*i+1], flteps, "ccopy(case 78) imag");
  58.      };
  59.    };
  60.   };
  61.  
  62.  
  63.   {
  64.    int N = 1;
  65.    double X[] = { 0.315, -0.324 };
  66.    int incX = 1;
  67.    double Y[] = { -0.312, -0.748 };
  68.    int incY = -1;
  69.    double expected[] = { 0.315, -0.324 };
  70.    cblas_zcopy(N, X, incX, Y, incY);
  71.    {
  72.      int i;
  73.      for (i = 0; i < 1; i++) {
  74.        gsl_test_rel(Y[2*i], expected[2*i], dbleps, "zcopy(case 79) real");
  75.        gsl_test_rel(Y[2*i+1], expected[2*i+1], dbleps, "zcopy(case 79) imag");
  76.      };
  77.    };
  78.   };
  79.  
  80.  
  81.   {
  82.    int N = 1;
  83.    float X[] = { 0.222f };
  84.    int incX = -1;
  85.    float Y[] = { 0.522f };
  86.    int incY = 1;
  87.    float expected[] = { 0.222f };
  88.    cblas_scopy(N, X, incX, Y, incY);
  89.    {
  90.      int i;
  91.      for (i = 0; i < 1; i++) {
  92.        gsl_test_rel(Y[i], expected[i], flteps, "scopy(case 80)");
  93.      }
  94.    };
  95.   };
  96.  
  97.  
  98.   {
  99.    int N = 1;
  100.    double X[] = { 0.021 };
  101.    int incX = -1;
  102.    double Y[] = { 0.898 };
  103.    int incY = 1;
  104.    double expected[] = { 0.021 };
  105.    cblas_dcopy(N, X, incX, Y, incY);
  106.    {
  107.      int i;
  108.      for (i = 0; i < 1; i++) {
  109.        gsl_test_rel(Y[i], expected[i], dbleps, "dcopy(case 81)");
  110.      }
  111.    };
  112.   };
  113.  
  114.  
  115.   {
  116.    int N = 1;
  117.    float X[] = { 0.376f, 0.229f };
  118.    int incX = -1;
  119.    float Y[] = { 0.143f, -0.955f };
  120.    int incY = 1;
  121.    float expected[] = { 0.376f, 0.229f };
  122.    cblas_ccopy(N, X, incX, Y, incY);
  123.    {
  124.      int i;
  125.      for (i = 0; i < 1; i++) {
  126.        gsl_test_rel(Y[2*i], expected[2*i], flteps, "ccopy(case 82) real");
  127.        gsl_test_rel(Y[2*i+1], expected[2*i+1], flteps, "ccopy(case 82) imag");
  128.      };
  129.    };
  130.   };
  131.  
  132.  
  133.   {
  134.    int N = 1;
  135.    double X[] = { -0.265, -0.84 };
  136.    int incX = -1;
  137.    double Y[] = { -0.156, 0.939 };
  138.    int incY = 1;
  139.    double expected[] = { -0.265, -0.84 };
  140.    cblas_zcopy(N, X, incX, Y, incY);
  141.    {
  142.      int i;
  143.      for (i = 0; i < 1; i++) {
  144.        gsl_test_rel(Y[2*i], expected[2*i], dbleps, "zcopy(case 83) real");
  145.        gsl_test_rel(Y[2*i+1], expected[2*i+1], dbleps, "zcopy(case 83) imag");
  146.      };
  147.    };
  148.   };
  149.  
  150.  
  151.   {
  152.    int N = 1;
  153.    float X[] = { 0.074f };
  154.    int incX = -1;
  155.    float Y[] = { -0.802f };
  156.    int incY = -1;
  157.    float expected[] = { 0.074f };
  158.    cblas_scopy(N, X, incX, Y, incY);
  159.    {
  160.      int i;
  161.      for (i = 0; i < 1; i++) {
  162.        gsl_test_rel(Y[i], expected[i], flteps, "scopy(case 84)");
  163.      }
  164.    };
  165.   };
  166.  
  167.  
  168.   {
  169.    int N = 1;
  170.    double X[] = { -0.374 };
  171.    int incX = -1;
  172.    double Y[] = { -0.161 };
  173.    int incY = -1;
  174.    double expected[] = { -0.374 };
  175.    cblas_dcopy(N, X, incX, Y, incY);
  176.    {
  177.      int i;
  178.      for (i = 0; i < 1; i++) {
  179.        gsl_test_rel(Y[i], expected[i], dbleps, "dcopy(case 85)");
  180.      }
  181.    };
  182.   };
  183.  
  184.  
  185.   {
  186.    int N = 1;
  187.    float X[] = { 0.084f, 0.778f };
  188.    int incX = -1;
  189.    float Y[] = { 0.31f, -0.797f };
  190.    int incY = -1;
  191.    float expected[] = { 0.084f, 0.778f };
  192.    cblas_ccopy(N, X, incX, Y, incY);
  193.    {
  194.      int i;
  195.      for (i = 0; i < 1; i++) {
  196.        gsl_test_rel(Y[2*i], expected[2*i], flteps, "ccopy(case 86) real");
  197.        gsl_test_rel(Y[2*i+1], expected[2*i+1], flteps, "ccopy(case 86) imag");
  198.      };
  199.    };
  200.   };
  201.  
  202.  
  203.   {
  204.    int N = 1;
  205.    double X[] = { 0.831, -0.282 };
  206.    int incX = -1;
  207.    double Y[] = { -0.62, 0.32 };
  208.    int incY = -1;
  209.    double expected[] = { 0.831, -0.282 };
  210.    cblas_zcopy(N, X, incX, Y, incY);
  211.    {
  212.      int i;
  213.      for (i = 0; i < 1; i++) {
  214.        gsl_test_rel(Y[2*i], expected[2*i], dbleps, "zcopy(case 87) real");
  215.        gsl_test_rel(Y[2*i+1], expected[2*i+1], dbleps, "zcopy(case 87) imag");
  216.      };
  217.    };
  218.   };
  219.  
  220.  
  221. }
  222.